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ABSTRACT 

Primordial magnetic fields (PMF) damp at scales smaller than the photon diffusion and free- 
streaming scale. This leads to heating of ordinary matter (electrons and baryons), which affects 
both the thermal and ionization history of our Universe. Here, we study the effect of heating 
due to ambipolar diffusion and decaying magnetic turbulence. We find that changes to the 
ionization history computed with recfast are significantly overestimated when compared 
with CosmoRec. The main physical reason for the difference is that the photoionization coeffi¬ 
cient has to be evaluated using the radiation temperature rather than the matter temperature. A 
good agreement with CosmoRec is found after changing this aspect. Using Planck 2013 data 
and considering only the effect of PMF-induced heating, we find an upper limit on the r.m.s. 
magnetic field amplitude of Bq ^ 11 nG (95% c.l.) for a stochastic background of PMF with 
a nearly scale-invariant power spectrum. We also discuss uncertainties related to the approxi¬ 
mations for the heating rates and differences with respect to previous studies. Our results are 
important for the derivation of constraints on the PMF power spectrum obtained from mea¬ 
surements of the cosmic microwave background anisotropies with full-mission Planck data. 
They may also change some of the calculations of PMF-induced effects on the primordial 
chemistry and 21cm signals. 
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1 INTRODUCTION 

The damping of primordial magnetic fields (PMF) heats electrons 
and baryons through dissipative effects ( [Jedamzik et al.|1998 [ |Sub-| 
|ramanian & Barrow|1998| l. This causes two interesting signals in 
the cosmic microwave background (CMB). One is due to the ef¬ 
fect on the CMB spectrum: the extra energy input from dissipat¬ 
ing PMF through the electrons leads to up scattering of CMB pho¬ 
tons, creating ay-distortion (see |Chluba & Sunyaev|20r^|Chlub^ 
|2014[ for recent overview on distortions) after recombination, with 
y-parameter up to y few x 10“^ i Jedamzik et al.| 2000 Sethi & 


|Subramanian||2005l [Kunze & Komatsu||2014| |2015|l. The second 

signal is seen as a change of the CMB anisotropies ( [Sethi & Sub-| 
|ramanian||2005[ [Kunze & Komatsu]|2014| [2bl5| l: the damping of 
PMF heats electrons above the CMB temperature. This reduces the 
effective recombination rate of the plasma, leading to a delay of re¬ 
combination and modifications of the Thomson visibility function. 


The changes of the CMB anisotropy power spectra caused by 
PMF-induced heating adds to the effects of PMF on the Einstein- 
Boltzmann system of cosmological perturbation^ The latter ef¬ 
fects have been subject of several investigations (see [Giovanninij 
|2006| for review) and were used to derive upper limits on the PMF 
amplitude smoothed on 1 Mpc scale of Bi Mpc S few x nG with 
pie-Planck [Paoletti & Finelli|201 ijjShaw & Lewis|2012[|Paoletti| 


& Finelli]|2013t and Planck ( [Planck Collaboration et ai.||2014c| 


2015[l data. In this paper, we discuss limits on the PMF ampli¬ 


tude using Planck 2013 data ( [Planck Collaboration et al.[2014a[b| (, 
only considering the effect of PMF-induced heating on the CMB 
anisotropies caused by changes in the ionization history. 

Previously, an approach similar to recfast jSeager et al.[ 
[2000|l was used to estimate the effects on the CMB energy spec¬ 
trum and CMB anisotropies (e.g., Kunze & Komatsu 2014[l. While 
for the standard cosmology ([Bennett et al.[[2003[ P anck Collab-[ 


[oration et al.|2014c^, recfast reproduces the calculations of de- 
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* In the standard treatment of magnetically induced cosmological pertur¬ 
bations, PMF contribute to energy density, pressure terms and and generate 
a Lorentz force on baryons. 
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tailed recombination codes like CosmoRec ( [Chluba & Thom^ 
201 1[( and HyRe c ( |Ali-Haimoud & Hirata|20Tl) very well ( jShaw 
& Chluba|201 1[ (, it was not developed for non-standard scenarios. 
We find that the effect of heating by PMF on the ionization his¬ 
tory is overestimated because the photoionization coefficient inside 
recfast is evaluated using the electron temperature, T^. Physi¬ 
cally, the heating does not alter the CMB blackbody significantly. 
Thus, the photoionization coefficient should be evaluated using the 
CMB blackbody temperature, as discussed in jChluba et al.jpQlO^ 
and also done in jSethi & Subramanian| ( [2()05| l. When comparing to 
CosmoRec, we find that after this modification the agreement be¬ 
comes very good. We also confirmed that collisional ionizations of 
both hydrogen and helium remain negligible until electron temper¬ 
atures Te - 10“* K at redshift z < fewx 10“ are reached. This is only 
achieved for significant PMF heating by strong fields, in excess of 
current limits on the magnetic field amplitude of a few nG. Thus 
the main effect that changes the ionization history, is a reduction of 
the recombination rate due to the higher electron temperature rather 
than extra ionizations. 

Our computations have two implications for the CMB. Firstly, 
the effect on the CMB power spectra from heating by PMF is re¬ 
duced significantly. This implies that expected limits on the spectral 
index and amplitude of PMF derived from CMB anisotropy mea¬ 
surements become weaker. Secondly, the departure of the electron 
temperature from the CMB temperature increases, because Comp¬ 
ton cooling is reduced. Still, we find that the difference in the y- 
distortion signal from heating by PMF is not as large, since reduc¬ 
tion of the free electron fraction and increase of the electron tem¬ 
perature more or less cancel each other (overall the same amount of 
energy is transferred to the CMB). The aspects discussed here may 
also be relevant to computations of the effect of PMF on the primor¬ 
dial chemistry (e.g., [Sc hleicher et al.||2 008a|b^ and 21c m signals 
(e.g., [Schleicher et al.[2009[[Sethi & Subramanian[2Q0f^ . Modified 
versions of CosmoRec and recfast-i-i-, which include heating by 
PMF will be made available at www. Chluba. de/CosmoRec 


2 COMPUTATIONS OF THE IONIZATION HI S TORY 


We consider a stochastic background of non-helical PMF charac¬ 
terized by 


(B,(k)S*(h)> = {Inf - h) Pij{k) PB(fc)/2, (1) 


where Pij(k) = dij - kjcj and P^{k) = Ak"^ determines the PMF 
power spectrum, with amplitude. A, and spectral index, «b- As 
a measure for the comoving integrated squared amplitude of the 
PMF, = Bq/A, we consider the treatment adopted by 
[Komatsu] ( [2014t : 


Kunze & 


/ 1 \ ^4. ^ 


2^1 


d(tr»+^ exp 

1 


^^O.b+3)/2 

2^f(^^)’ 


where a Gaussian filter is used and the damping scale is given by 
[Kunze & Komatsu] [2014[ l: 

ku X 286.91 (nG/Bo)Mpc-‘. (3) 


We use the CGS conventions, i.e. (ps) = B^ISn is the average co¬ 
moving magnetic field energy density. 

To include the effect of dissipation from PMF, we follow the 
procedure of [Sethi & Subramanian[ ( |2005) and [Kunze & Komatsu] 
( [2014[ henceforth, KK14). We implemented both the heating (F = 
dE I dl = released energy per volume and second) by ambipolar 


diffusion and decaying magnetic turbulence (see Appendix [A| for 
additional details). The PMF heating has to be added to the electron 
temperature equation with 


dTe 

df 


= -2Hff + 


ScTTNe Py 
'im^cNiot 


(Ty - Te) -t 


F 

(3l2)kN,,,' 


(4) 


Here, H{z) denotes Hubble rate, Atot = Ah(1 + /hc + Ae) the num¬ 
ber density of all ordinary matter particles that share the thermal 
energy, beginning tightly coupled by Coulomb interactions; Ah is 
the number density of hydrogen nuclei, /ne ~ Fp/4(1 - Tp) 0.079 
for helium mass fraction Fp = 0.24; = NffNu denotes the free 

electron fraction and Py = affP* x 0.26 eV(l -i- zf the CMB energy 
density. The first term in Eq. 0 describes the adiabatic cooling 
of matter due to the Hubble expansion, while the second term is 
caused by Compton cooling and heating. The last term accounts 
for the PMF heating. Notice that the last two terms in Eq. 0 differ 
slightly from those presented in earlier works (e.g., [Sethi & Sub-| 
[ramanian poos] ). One reason is that the heat capacity contribution 
from helium was neglected so that for the Compton cooling term 
Mot ~ Ah( 1 -1-Ae). Secondly, for the PMF-induced heating term, the 
thermal energy was distributed only among the hydrogen atoms, 
Atot X! Ah = Ae/Ae, although even without helium, the free elec¬ 
trons contribute. However, we find that this only changes the free 
electron number by ANffN^ x 10% - 20% in the freeze-out tail. 

We modified both recfast-i-+ (which by default is meant to 
reproduce the original version of recfast) and CosmoRec to in¬ 
clude the effects of heating by magnetic fields. For recfast+-i-, we 
can separately adjust the computation of the photoionization coef¬ 
ficients for hydrogen and helium. In the default setting, they are 
evaluated using the matter temperature, T = T^, as in recfast. 
From the physical point of view, the photoionization coefficient, 
jS/c, of an atomic level i is a function of both photon and elec¬ 
tron temperatures, Ty and T^, respectively. The dependence on 
the electron temperature enters through Doppler boosts. Even for 
high electron temperatures, this correction can be neglected, so 
that one has/3,c = fiic(Ty,ff) x ff^iTy) x 4trf ^^^<T,e(v) dv, 
where By(Ty) is the CMB blackbody intensity (e.g.,|Seager et al.j 
[2000[ (. Clearly, without significant distortions of the CMB radia¬ 
tion field, this expression directly confirms that the photoionization 
rate, B,c = A,-jS/cfTy), where A, is the population of the level i of the 
atom, depends only on the photon temperature. Thus, the effective 
(case-B) photoionization rate also only depends on the photon tem¬ 
perature, a modification that causes a big difference for the effect 
of heating by PME, as we show below. 

In contrast to this, the photo-recombination coefficient, Od, to 
atomic level i mainly depends on the electron temperature, with a 
smaller correction due to stimulated recombinations in the ambi¬ 
ent CMB blackbody radiation field. This implies, ad = UdiT^, Ty), 
which for the recfast treatment is set to ad ~ ad(T,., Ty = Ty.). 
For the detailed recombination calculations, this approximation 
becomes inaccurate for highly excited levels (e.g., [Chluba et al.[ 
[2007|, changing the freeze-out tail of recombination at the per¬ 
cent level jChluba et al.|2007l [Grin & Hirata|2010[[0iluba et ah] 
[2010[ l. In CosmoRec and HyRec, the full temperature dependence 
of the photo-recombination coefficient is taken into account us¬ 
ing an effective multi-level atom method (jAli-Haimoud & Hirataj 
[2010[ (. When including heating by PMF, the recfast treatment 
thus slightly overestimates the photo-recombination rate to each 
level, since for Ty « T^. stimulated recombinations are overesti¬ 
mated when assuming Ty - T^. However, the difference is much 
less importanf fhan the error caused by evaluating the photoioniza¬ 
tion coefficient for T = T,.. 
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Magnetic heating 3 


2.1 Collisional ionization 


The exponential dependence on the ionization potential suppresses 
the effect of collisional ionization from the ground state, so that 
in the standard computation they can be neglected ^Chluba et al.| 
|2007| >. Since at low redshifts (z < 800) the electron temperature 
can be pushed quite significantly above the CMB photon tempera¬ 
ture by heating processes (lower panels in Fig.[^an d|2), it is impor¬ 
tant to check if collisional ionizations by electron impacpjbecome 
efficient again. Using the fits of |Bell et d^ ( |1983[ l 


dA'is.m 

df 

dMs.Hel 


« -5.85 X lO-'^r^'^e-^H/^cm^ s^* Ms.niA'e 
« -2.02 X lO^^’r^'^e^^H'/^cm^ s^‘ Al|s,Hei A^e 


with Th ~ 1.58 X 10 ^ K, Tne ~ 2.85 x 10 ^ K and T^ = 7 ’/ 10 ‘* K, 
where T = T^, we confirm that this effect can usually be neglected. 
We nevertheless add these rates to the calculation whenever heat¬ 
ing by PMF is activated and for very large heating (pushing the 
electron temperature up to Te - 10** K) they do become important 
in limiting the maximal electron temperature. We also included the 
cooling of electrons by the collisional ionization heating to ensure 
the correct thermal balance. 


2.2 Decaying magnetic turbulence 

Using recfast++ with default setting, we are able to reproduce 
the central panel in Fig. 10 of KK14 for decaying magnetic turbu¬ 
lence. One example, for Bq = 3 nG and = -2.9 is shown in 
Fig. □ We compare the standard recombination history (no extra 
heating) with three cases obtained from recfast+-t and the full 
computation of CosmoRec. The effect of reionization at z < 10 was 
not included (see |Kunze & Komatsu|[20T5| for some discussion), 
as it does not affect our main discussion. The first agrees well with 
the result of KK14, with a large change in the freeze-out tail of 
the recombination history being found (dotted line). Modifying the 
evaluation of the hydrogen photoionization rate to T = Ty gives a 
smaller change (dash-dotted line). Also changing the evaluation of 
the helium photoionization rate finally gives the dashed line, with 
a 5 times smaller effect on the freeze-out tail. Using the standard 
recfast approach, the photoionization rate is thus overestimated 
so that even helium is partially reionized. We find that after chang¬ 
ing the evaluation of the photoionization rates lo T = Ty the result 
obtained with recfast++ agrees to within =: 10% with the detailed 
treatment of CosmoRec (solid/red line). This case is also fairly close 
to the result for m = 2(nB + 3)/(nB -t- 5) = 0.1(= hb - -2.9) shown 
in Fig. 4 of |Sethi & Subramanian| ( (2()05| >. The remaining difference 
to CosmoRec is caused by stimulated recombination effects that are 
not captured correctly with a recfast treatment. 

Our computations show that the smaller effect on the free elec¬ 
tron fraction allows the electron temperature to rise higher above 
the photon temperature than with the default recfast treatment 
(see Fig. [TJ. This is because for a lower free electron fraction, 
Compton-cooling becomes less efficient. We find that in terms of 
the Compton-y parameter, these two effects practically cancel each 
other, leaving a difference at the level of < 5%. For instance, com¬ 
puting the y-parameter, y = / ttyAeC dt, for Bq - 3nG and 

ns = —2.9 using the default recfast-i-i- result we obtain y 
1.0 X 10“^(So/3 nG)^, while when evaluating the photoionization 


^ Protons are heavier and thus slower, so that their effect is much smaller. 




Figure 1. Effect of heating by decaying magnetic turbulence on the ion¬ 
ization history (upper panel) and electron temperature (lower panel) for 
Bq = 3 nG and ns = -2.9. 


rates correctly we have y == 9.7 x 10“*, corresponding to a 4% ef¬ 
fect. For larger spectral index, the difference becomes even smaller. 
For Bq = 3nG and ns = 0, we find y == 5.4 x 10“^(Bo/3 nG)^ 
with a difference < 1% in the two treatments. The reason is that 
for larger spectral index, most of the effect arises from higher red- 
shifts (z - 10*), which are less sensitive to the evaluation of the 
photoionization rates since Compton-cooling is still extremely effi¬ 
cient, forcing T^=^Ty. 

In the treatment of the heating by decaying magnetic turbu¬ 
lence, we switch the effect on rather abruptly (Az/z - 5%) around 
Zi — 1088 following previous approaches i Sethi & Subramanian| 


|2005[|Schleicher et al.|2008b[|Kunze & Komatsu|2014| . Although 

the effect of heating by decaying magnetic turbulence is not as 
visible at early times (see Fig. [T](, this approximation adds un¬ 
certainty to the predictions of the CMB anisotropies since small 
effects close to the maximum of the Thomson visibility function 
can have a larger effect than similar changes in the freeze-out tail 
(e.g., |Rubino-Martm et al.|2dO^|Farhang et al.|2012^ . For detailed 
CMB constraints, this approximation should be improved, includ¬ 
ing more detailed consideration of the time-dependence of the heat¬ 
ing rate at z > z, - 1088. For example, when changing from very 
abrupt to more smooth transition between no heating and heating 
at z - Zi, we find that the numerical result for the TT power spec¬ 
trum at large scales {{ < 200) is affected noticeably. However, in 
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Figure 2. Effect of ambipolar diffusion on the ionization history (upper 
panel) and electron temperature (lower panel) for Bq = 3 nG and = -2.9. 
For comparison, we also show the case obtained with settings similar to 
[Sethi & SubramanianH20()5) . 


the present paper we shall follow the previous approach, addressing 
order-of-magnitude questions only. 


2.3 Ambipolar diffusion 

For heating by ambipolar diffusion, we also find a reduction of the 
effects when modifying the evaluation of the photoionization rates. 
This is illustrated in Fig.[^ again for Bq = 3 nG and ng = -2.9. The 
effect of reionization was not modeled. Evaluating the photoion¬ 
ization rates correctly reduces the effect on the freeze-out tail by 
more than one order of magnitude. Also, the low-redshift electron 
temperature is underestimated by roughly one order of magnitude, 
with a plateau rising to Te - lO'^ K. Comparing with Figures 1 and 
2 of [Sethi & Subramanian| ( [2005) , we find good agreement with 
our computation when setting x p\ll\, which is equivalent 
to setting /l = 1 (see lines labeled ‘SS2005’ in our Fig.|^ for the 
average Lorentz force caused by the PMF. 

With the expressions given in Appendix]^ we are also able to 
reproduce Fig. 10 of KK14 for the ambipolar diffusion case; how¬ 
ever, we had to multiply our heating rate by a factor of l/(8;r)^, 
albeit being based on the same expressions. We confirmed the or¬ 
der of magnitude of the heating rates for ng = -2.9 using an alter¬ 
native evaluation based on the analytic expressions of [Finelli et al.j 


( [2008^ and [Paoletti et al.[f2009^ , finding that the importance of the 
process was indeed underestimated. This changes the y-parameter 
contribution caused by ambipolar diffusion. Including only heating 
due to ambipolar diffusion for Bq = 3 nG and ng = -2.9 we find 
y X 3.1 X 10“*'(5o/3 nG)^ for the heating rate of KK14, while here 
we findy =: 6.6x 10“®(6o/3 nG)^. Relative to the contribution from 
decaying magnetic turbulence this is a small correction, well be¬ 
low the precision of the approximations made in the computation. 
However, for larger spectral index, the effect becomes more impor¬ 
tant. For Bq - 3 nG and ng = 0 we find y =: 7.5 x 10“**’(Bo/3 nG)^ 
using the evaluation of KK14 but y x Fix 10“’(Bo/3 nG)^ with 
our treatment, making this contribution comparable to the one from 
heating by decaying magnetic turbulence. 

For ambipolar diffusion, the numerical treatment mainly de¬ 
pends on the evaluation of the average of the Lorentz force squared 
and the distribution used to characterize the stochastic background 
of PMF as small scales. [Sethi & SubramanianH2005[ l and [Schl^ 
[icher et ^ ([2008b), adopted an order of mag nitude estimate , i.e.. 


Kunze 


1 2011 1 were 


X while in KK14, expressions from 
used. This introduces a strong dependence of the heafing by am¬ 
bipolar diffusion on ng, reducing the effect significantly as the spec¬ 
tral index approaches ng == -3. Using a sharp cut-off instead to 
approximate the effect of damping at small-scales allows obtaining 
exact analytic expressions for the energy-momentum tensor corre¬ 
lators ( [Finelli et al.|2008[ [Paoletti et al.|2009T l. However, here we 
adopt the approximations of[Kunze]j2011^ to illustrate the effects. 


3 EFFECTS ON THE CMB ANISOTROPIES 

The changes to the ionization history introduced by heating from 
magnetic field inevitably affects the CMB temperature and polar¬ 


ization anisotropies jSethi & Subramanian[[2005] [Kunze & Ko-[ 


matsu|2014| 

2015J. We can estimate the importance of this effect 

using camb 

Lewis et al. [[20001 with our modified recombination 


codes. We use the Planck 2013 cosmology ([Planck Collaboration[ 
[et al.[[^14c[ l with reionization optical depth t = 0.09. Changes 
of the free electron fraction around decoupling z == 1100 usually 
weigh more than modifications at late times in the freeze-out tail 
(e.g. [Rubino-Martm et ar][2008[ [Farhang et al.[[201^ . Thus, al¬ 
though not as visible in Fig. and [^ especially at small scales, 
a large part of the effect on the CMB power spectra arises from 
modifications z - 1100, causing additional diffusion damping 
and shifts in the positions of the acoustic peaks. 

For ambipolar diffusion, the effect on the ionization history 
around z x 1100 is much smaller than for decaying magnetic tur¬ 
bulence. The former mainly affects the freeze-out tail and thus the 
optical depth to the last scattering surface, t. This leads to extra 
e^^’’ damping of the CMB anisotropies at small scales, an ef¬ 
fect that is partially degenerate with the curvature power spectrum 
amplitude, A^, and its spectral index n^. Extra polarization at large 
scales is generated by re-scattering events, an effect to which CMB 
polarization data is sensitive (see also [Kunze & Komatsu[[20T5) (. 
These effects are very similar to changes to the CMB power spectra 
caused by dark matter annihilatio n ([Chen & Kamionkowski|2004[ 
[Padmanabhan & Finkbeiner[2005[[Zhang et al.[2006^ , so that some 
degeneracy with this process is expected, especially when includ¬ 
ing the effect of clumping at late times ( [Hiitsi et al.|2009^ , which 
can boost the free electron fraction in the freeze-out tail in a similar 
manner. For decaying magnetic turbulence, the Thomson visibil¬ 
ity function is affected close to its maximum, so that changes in 
the positions of the acoustic peaks are found, which can be tightly 
constrained using CMB data. 
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Figure 3. Effect of PMF heating on the TT and EE power spectra for two 
choices of the PMF power spectrum parameters. The standard Planck 2013 
cosmology with r = 0.09 was used. See Sect. I^for discussion. 

In Fig. 1^ we show the separate contributions from decaying 
magnetic turbulence and ambipolar diffusion to the changes in the 
TT and EE power spectra. For Bq = 3 nG and = -2.9, the 
effect of ambipolar diffusion is small and the dominant effect is 
caused by decaying magnetic turbulence, which introduces clear 
shifts in the positions of the acoustic peaks. Setting = 3.0, 
we see that ambipolar diffusion does add a significant correction. 
At -6%(So/3nG)^, to the Thomson optical depth. At very 
low t, it also introduces features into the power spectra, which 
help breaking the degeneracies mentioned above. The shifts in the 
peak positions due to decaying magnetic turbulence also increase 
strongly for this case. Notice that the heating rates for both decay¬ 
ing magnetic turbulence and ambipolar diffusion scale as F oc Bg, 
so that the changes in the CMB power spectra strongly decrease 
with Bq. The aforementioned effects can be constrained with cur¬ 
rent data, and one expects decaying magnetic turbulence to drive 
the limits, at least for quasi-scale invariant PMF power spectra. 



0.0 I-^^ -- ' - 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 


Bq [nG] 

Figure 4. Posterior distribution for the PMF strength, Bq, with = -2.9 
and different combinations of heating by decaying magnetic turbulence 
(MHD) and ambipolar diffusion. 

the Einstein-Boltzmann system for cosmological perturbations, the 
Planck 2013 95% c.l. upper limit on the magnetic field strength, 
smoothed over 1 Mpc length, if]Si Mpc ^4.1 nG, obtained by vary¬ 
ing riB in the interval [-2.9,3]. 

In Fig. [^ we show the posterior distribution for Bq at fixed 
= -2.9. We compare three cases, including the heating caused 
by decaying magnetic turbulence, ambipolar diffusion and the com¬ 
bination of both. The 95% upper limits on the magnetic field 
strength are < l.lnG, < 1.5nGandSo :S 1.1 nG, 

respectively. As anticipated earlier, decaying magnetic turbulence 
strongly drives the constraint and shape of the posterior distribu¬ 
tion of Bq, with ambipolar diffusion leading to a correction only. It 
is also clear that the limits 5“*' and Bq are so comparable 

mainly because the posteriors have strong non-Gaussian tails. 

Evaluating the photoionization rates using T = T^, we expect 
an upper limit of Bq < 0.5 nG (95% c.l.) when including both am¬ 
bipolar diffusion and decaying magnetic turbulence. This is about 
- Vs - 2.2 times tighter than the limit quoted above, simply be¬ 
cause the effect on the ionization history is overestimated. This il¬ 
lustrates how important the modification to the recfast treatment 
discussed here is. 


4 CONCLUSIONS 

We investigated the effect of heating due to ambipolar dijfusion 
and decaying magnetic turbulence on the thermal and ionization 
history of our Universe. We find that changes in the ionization his¬ 
tory, computed with an approach similar to recfast, are signif- 


3.1 Constraints from Planck 2013 data 

In this section, we discuss constraints on the PMF power spectrum 
using Planck 2013 data (Planck Collaboration et al.|20I4c^ . We ex¬ 
plicitly include only the effect of PMF-induced heating on the CMB 
power spectra. Taking into account only the PMF contributions to 


^ In |Planck Collaboration et ^(2014^ the constraint on the amplitude of 
PMF is quoted in terms of Bi Mpc, i e., the amplitude smoothed over 1 Mpc 


Shaw & Lewis 

2012 

have therefore 

DO - 


tiB = -2.9. 


. Note that the relation B^ = (Arod/ 

1.3 Bi Mpc for kd given by Eq. |^, A = I Mpc and 
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icantly overestimated when compared with CosmoRec. However, 
after evaluating the photoionization rates at the photon tempera¬ 
ture, T = Ty, the results agree to within =: 10% with the more de¬ 
tailed treatment. The remaining difference to CosmoRec is mainly 
because of stimulated recombination effects. 

For the Compton y-parameter computed in different treat¬ 
ments of the problem, we find only small differences at the level 
of < 5% for decaying magnetic turbulence. This is because the 
reduction of the effect on the free electron fraction is roughly com¬ 
pensated by the increase in the electron temperature. However, our 
computations do show that the y-parameter contribution caused by 
ambipolar diffusion was underestimated. For nearly scale-invariant 
PMF power spectrum, ambipolar diffusion still causes only a small 
correction relative to the y-parameter contribution from decaying 
magnetic turbulence; however, for «b - 0 the two contributions 
become comparable in order of magnitude. 

Using Planck 2013 data and only including the PMF heat¬ 
ing effect, we find an upper limit on the magnetic field strength 
of Ba < 1.1 nG (95% c.l.) for a PMF power spectrum with spectral 
index wb = -2.9. As shown in |Planck Collaboration et aLl ( |2015| l, 
the heating effect considered here leads to a tighter constraint than 
the one derived by considering only the direct effects of PMF on the 
cosmological perturbations. The improvement is approximatively a 
factor of 3 for Planck 2015 data (Planck Collaboration et al.|2015t 
and hb = -2.9. However, we expect uncertainties in the model¬ 
ing of the heating by decaying magnetic turbulence to affect the 
results for quasi-scale invariant PMF power spectra, while for blue 
PMF power spectra, details in the modeling of ambipolar diffusion 
becomes important. Given that the effects of PMF heating studied 
here are so competitive with those of PMF on the fluid perturba¬ 
tions, this problem deserves more careful consideration, in partic¬ 
ular with respect to future improved measurements of polarization 
on large angular scales by Planck, which could lead to a better es¬ 
timate of T. 

An additional uncertainty is caused by the way the reioniza¬ 
tion epoch is added at z < 10. Currently, we simply use the camb 
default prescription that ensures a smooth transition at the start 
of reionization. However, for strong PMF, the pre-reionization at 
z > 10 can be significant (see Fig. [T] and |^, so that priors on the 
total optical depth need to be considered more carefully. This is ex¬ 
pected to be important only for very blue PMF power spectra or 
when the PMF amplitude becomes very large. 

For quasi-scale invariant PMF power spectra, our analysis sug¬ 
gests that PMF heating can contribute no more than y < 1.1 x 10“* 
(95% c.l.) to the average Compton y distortion. Although measure¬ 
ments with a PIXIE-like experiment (Kogut etal.|2011^ could reach 
this sensitivity, a much larger distortion (y =; 10“^ - 10“*) is cre¬ 
ated just from the reionization and structure formation process (Hu| 
|et al.|199^|Cen & Ostriker|1999[[^fregier & et al.|200^ . Thus, 
it will be very difficult to use future spectral distortions measure¬ 
ments to constrain the presence of PMF in the early Universe in a 
model-independent way. 

Finally, we confirmed that evaluating the photoionization rates 
at T = Ty in recfast++ does not affect the ionization history by 
more than AAe/Ae - 0.3% at z - 780 for the standard cosmol¬ 
ogy and thus has no significant effect on the analysis of current and 
upcoming CMB data. To improve the consistency of the recfast 
treatment one could thus change this convention and then recali¬ 
brate the fudge-functions without additional changes. 

Our results are important for the derivation of constraints on 
the PMF power spectrum obtained from measurements of the cos¬ 
mic microwave background anisotropies with Planck full-mission 


data (see |Planck Collaboration et al.|2015') ). They may also be rele¬ 
vant to computations of the effect of PMF on the primordial chem¬ 
istry (e.g., Schleicher et al.|2008a|b^ and 21c m signals (e.g., |Schi^ 
|icher et ah 2009||Sethi & Subramanian|2009t . In all cases, it will be 
important to improve the description of the PMF heating rate, since 
current approximations introduce noticeable uncertainty. This will 
be left for future work. 
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Magnetic heating 1 


APPENDIX A: HEATING EROM DECAYING MAGNETIC 
TURBULENCE AND AMBIPOLAR DIEEUSION 

To describe the heating caused by magnetic fields we follow the proce¬ 
dure of|Sethi & Subramanian| j2005) , with some specific parameterizations 
given in |Schleicher et al.H2QQ8a) '^d KK14. For more discussion about the 
physics of the problem, we refer to these references. 


A1 Decaying magnetic turbulence 


The heating rate caused by decaying magnetic turbulence can be approxi¬ 
mated as ISethi & Subramanian|2QQ5) 

fln(l + f)r 

-ym.z)p^{z\ (Al) 




with m = 2(nB + + 5), tj/td ~ 14.8(So/nG)“VA:d/Mpc“^)“^ damp¬ 

ing scale kd « 286.91 (Bo/nG)“^Mpc“\ and magnetic field energy den¬ 
sity pb(G = ^0*^1 +z)^l{^n) ^ 9.5 X \Q~^(Bq/ nG)^ py{z), where py{z) « 
0.26 eV cm“^(l + z)‘^ is the CMB energy density. In the approximation, the 
heating switches on abruptly at redshift Zi = 1088, however, in our compu¬ 
tation we switch the heating on more smoothly to avoid numerical issues. A 
refined physical model for the heating is required to improve this treatment. 


A2 Ambipolar diffusion 


To capture the effect of heating by ambipolar diffusion we use ^Sethi &| 
|Subramanian|2QQ31|Schleicher et al.|2008ai 


^ (l-Ap) (l(VxB)xBp) 

rXpPb 


(A2) 


where = ^|(V X B) X /(4;r)^ denotes the average square of the 
Lorentz-force andpb = the baryon mass density with baryon number 

density Ab. Neglecting corrections from helium, we only need the free pro¬ 
ton fraction, = Ap/An. to describe the coupling between the ionized and 
neutral component. The coupling coefficient is given by y = {cru)^^+ /2mH 
with (try)////+ « 6.49 X 10“^*^(7’/K)^-^^^cm^ s“^ 

For-3 < «B < 5, the integral for the Lorentz force, Eq. (3.5) of KK14, 
is well-approximated by 


(|(V X B) X Bp) « /i(„B + 3) = 16;rV|(z)^d"(z) A("b + 3) (A3) 


ftix) = 0.8313[1 - 1.020 X 


(A4) 


Here, B = Bo(l + and 1^ = a/ki. 
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